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In this work, the adsorption-desorption kinetics in the framework of the lattice gas model are 
analyzed. The transition probabilities are written as an expansion of the occupation configurations. 
Due to that the principle of detail balance determine half of the adsorption, At and desorption, 
Di coefficients, different functional relations between them are proposed. Introducing additional 
constrains, it is demonstrated that when those coefficients are linearly related through a parameter 
7, there are values of lateral interaction, V, that lead to anomalous behavior in the adsorption 
isotherms, sticking coefficients and thermal desorption spectra. Diagrams for the allowed values of 
V and 7 are also shown. Alternatively, a non-linear relation among those coefficients is introduced. 
In such way the equilibrium and non equilibrium observables do not present any anomalous or 
inconsistent behavior. 

PACS numbers: 02.50.-r, 68.43.Mn, 68.43.Nr, 47.11. Qr, 68.43.Vx 



I. INTRODUCTION 

The importance of surface and interface in a wide va- 
riety of fields has inspired an enormous interest in their 
structure and dynamics. The study of the kinetics and 
dynamics of the surface processes is of fundamental in- 
terest in the understanding of heterogeneous catalysis 
and other processes ta king place i n g as-solid interfaces 
0, i 1 i, & i, 0, ! , I, MM El M • The kinetics 
is fundamentally determined by the energy transfer [lj|. 

The rate of change in the adsorbate can be written 
as a difference between adsorption and desorption terms. 
This can be made in a heuristic way or in the frame- 
work of a rigorous approach using the non-equilibrium 
thermodynamic theory 0, [HI, [III • 

If the adsorbate does not remain in a quasi-equilibrium 
state during desorption, then the description of the sys- 
tem through macroscopic variables is not enough and it 
should be carried out based on the non-equilibrium statis- 
tical mechanics that involves time-dependent distribution 
functions. The Kinetic Lattice Gas Model (KLGM) is an 
example of such theory, and it was set up in close analogy 
to the time dependent Ising model for magnetic systems, 
which was originally introduced by Glauber [UHl]. In 
its simplest form the KLGM is restricted to the submono- 
layer regime and to the gas-solid system where the sur- 
face structure and the adsorption sites do not change 
with the coverage [l9j . However, further generalizations 
of the lattice-gas model have been made in the past, for 
example: the multilayer lattice gas model was used for 
analyzing and simulating desorption data for Li and Cs 
on Ru(0001) [13, HH; the influence of the precursor state 
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in the adsorption-desorption kinetics [161 ] and other re- 
lated problems [22|; the adsorption-desorption kinetics 
with multiple site occupations [lH, [24|, [25| , etc. 

On the other hand, it is well known that different mi- 
croscopic dynamics can yield different equilibrium paths 
and equilibrium fluctuations [2t| (cluster versus local 
MC algorithms being the most extreme example [271 ]) 
and even noticeable differences in the steady-state mi- 
crostructure [28l . [29j . Nevertheless, the general expec- 
tation is that, if no additional parameters (such as an 
activation barrier or a diffusion rate) are introduced into 
the physical model, observables are only affected quanti- 
tatively. 

Recent studies indicate that different stochastic dy- 
namics, even when they have the same conserved quan- 
tities and satisfy detailed balance, lead to important dif- 
ferences in the nanostructure of field-driven interfaces 

[M HE H2, m, [11] . 

Before defining the transition probabilities, an impor- 
tant distinction must be made. That is, between mod- 
els with hard dynamics [13, H [H, [H, HE H, HI, in 
which the single-site transition rates cannot be factor- 
ized into one term that depends only on the interaction 
energies, in contrast with those models with soft dynam- 
ics (all, M, m El, E3, S 0, EE m , for which this 

factorization is possible. 

In this context, Kang et al. (44| have analyzed differ- 
ent choices for the transition probabilities in MC simu- 
lations for studying the growth exponent in the domains 
growth. They showed that the choice of transition prob- 
abilities affects directly the dynamic quantities. Particu- 
larly, they found this dependence in the growth exponent 
in the ordered domains out of equilibrium [47| . In the 
same way, Rikvold et al. have shown that the intrinsic 
interface width and properties in field-driven solid-on- 
solid interface studies depend on the choice of dynamics. 
They have found that, in the framework of the soft dy- 



2 



namics, all dependence on the field canceled due to the 
principle of detailed balance [U [H, [Hj] . On the other 
hand in hard dynamics the intrinsic interface width and 
properties, such as the propagation velocity, are strongly 
affected by the field [p. 

Despite the fact that, detailed microscopic mechanisms 
of the surface processes are usually not known, in the 
KLGM the transition probabilities can be written in 
terms of the o ccup ation configurations of all neighbor- 
ing sites 0, [H, HI, [10, [ML 123, [H, Si] . In this point, 
the principle of detailed balance imposes a set of restric- 
tions on the coefficients of adsorption Ai, desorption Di, 
diffusion d, etc. as is discussed by H. J. Kreuzer and 
co-workers [IH, [ljl Ojl, [22[ , where the authors have intro- 
duced different kinetics according to the relations among 
those coefficients. When adsorption and desorption are 
the only processes taken into account, they have assumed 
linear relation Ai = "fDi (where 7 is a proportional- 
ity parameter). In such way the authors have obtained 
the adsorption-desorption kinetics calculating the stick- 
ing coefficients for different cases, including Langmuir ki- 
netics, sticking on random adsorbate and the influence of 
intrinsic and extrinsic precursors ,16]. Moreover, very re- 
cently, S. H. Payne and H. J. Kreuzer have discussed the 
one-dimensional diffusion with different lateral interac- 
tion, using the same methodology |4§ |. 

Although those coefficients satisfy the principle of 
detailed balance, the relation between them depends 
strongly on the lateral interaction energy and must be 
carefully analyzed according to the adsorption and des- 
orption processes, otherwise the results could be not gen- 
eral leading to anomalous behaviors in the kinetics. 

The aim of the present paper is to demonstrate that, 
even for the simplest one-dimensional linear re- 

lation among the adsorption and desorption coefficients 
gives some inconsistencies. With this purpose, the ad- 
sorption desorption kinetics are analyzed for a ID KLGM 
with nearest neighbor lateral interaction, V. 

The analytical treatment have been done in a one- 
dimensional system, mainly due to two reasons. The 
first one, is the availability of the exact solution for the 
coverage and two-sites correlation function (which is not 
possible in higher dimensions). The second one, is the 
presence of phase transitions in two and higher dimen- 
sions. This could mask the possible anomalous behaviors 
in the observables, which is precisely the objective of the 
present investigation. 

The outline of the paper is as follows: In Section II, 
the KLGM is set up, introducing the master equation 
and writing the transition probabilities in terms of the 
occupation configurations of all neighboring sites. A set 
of restrictions on the adsorption and desorption coeffi- 
cients are obtained according to the principle of detailed 
balance. In order to fix the unknown coefficients in the 
transition probabilities, two different linear relations are 
proposed. The equilibrium and non equilibrium proper- 
ties of the system are analyzed solving the rate equations 
for the coverage and higher correlations. Additional in- 



equalities among the adsorption and desorption parame- 
ters are found, determining forbidden values for the pa- 
rameters 7 and V, for which the adsorption isotherms, 
sticking coefficients and thermal desorption spectra are 
ill-behaved. Diagrams with forbidden regions as function 
of lateral interactions V and the parameter 7 are showed. 
In order to solve these inconsistencies, in Section IV, a 
functional relation between the adsorption and desorp- 
tion parameters is introduced. In such way the observ- 
ables present a correct behavior without inconsistencies. 
Finally, in Section V our conclusions are presented. 



II. THE KINETIC LATTICE GAS MODEL AND 
THE MASTER EQUATION 

To set up the KLGM one restricts the analysis to a gas- 
solid system in which all relevant processes, like diffusion, 
adsorption, desorption, reactions, etc., are Markovian. 
One assumes that the system can be divided into cells, 
labeled i, for which one introduces microscopic variables 
iii = 1 or depending on whether cell i is occupied by 
an adsorbed gas particles or not. The connection with 
magnetic systems is made by a transformation to spin 
variables <Xj = 2m — 1. To introduce the dynamic of the 
system one writes down a model Hamiltonian 

H = E a ^2m + -V^2n i n j + ... (1) 

' (ij) 

Here E s is a single particle energy, V is the two particle 
interaction between nearest neighbors (ij). Interactions 
between next nearest neighbors etc., and many particle 
interactions can be easily added to eq. (1). 

As long as the number of particles in the adsorbate 
does not changed, which is the case for systems in equilib- 
rium or diffusion studies, the first term in (1) is constant 
and can be dropped from further consideration. How- 
ever, if the objective is the study adsorption-desorption 
kinetics, the number of particles in the adsorbate changes 
as a function of time and a proper identification of E s is 
mandatory. Arguing that the lattice gas Hamiltonian 
should give the same Helmholtz free energy as a micro- 
scopic Hamiltonian (for noninteracting particles) one can 
show that the proper identification is given by [50( , 

E s = -V - \ Hq 3 q lnt ) - \ [ln(/3A 3 P) - ln(Z inl )] (2) 

P P 

where (3 — 1/ksT; ks and T are the Boltzmann con- 
stant and the absolute temperature, respectively. Vq is 
the (positive) depth of the surface potential, q 3 is the 
single particle partition function of an adsorbed parti- 
cle, qi n t is the internal partition function for (frustrated) 
vibrations and rotations of the adsorbed molecule; P 
is the pressure in the gas phase above the surface and 
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A = h/ \J2-KmkBT is the thermal wavelength of the ad- 
particle with mass m, which partition function for (free) 
vibrations and rotations is Zi nt . 

One introduces a function P(n, t) which gives the 
probability that a given microscopic configuration n = 
(rii, ri2, «at) is realized at time t, where TV is the total 
number of adsorption sites on the surface. It satisfies a 
master equation, 

dP(n; t) = J2 n)P(n; t) - W(n] n)P(n; t)] (3) 

n 

where W(n; n) is the transition probability that the 
microstate n changes into n' per unit time. It satisfies 
detailed balance 

W(n; n)Po(n) = W(n; n)P {n) (4) 

where 

Po(n) = Z- 1 e X p{-pH(n)) (5) 

is the equilibrium probability and Z is the partition 
function given by 

Z = Y / ^ fm{n) (6) 

In principle, W(ri; n) must be calculated from a Hamil- 
tonian that includes, in addition to (1), coupling terms to 
the gas phase and the solid that mediate mass and energy 
exchange. However, depending of the system, different 
expressions for transition probabilities can be proposed. 
In transition dynamic approximation [4ll |45| , transition 
rates cannot be factorized into one part that depends 
only on the interaction energy and another that depends 
only on the field energy [34j (hard dynamic (35l|). 

I 

dNN 

dt 
dNEN 

dt 

dNNN 
Jt 



The principle of detailed balance imposes a set of re- 
strictions on the coefficients A and Di [22 |. which are: 



Usually, one follows the procedure introduced by 
Glauber and guesses an appropriate form for W(n;n). 
One further assumes that the duration of an individual 
transition, e.g., hopping to a neighboring site, is much 
shorter than the residence time in the initial state. In 
this situation there will be only one transition at any 
given time and the total transition probability as a sum 
of individual terms can be written. 

In order to analyze the simplest cases, let us consider 
the one-dimensional lattice gas with nearest-neighbors in- 
teractions where only direct adsorption and desorption 
processes are taken into account (no other processes are 
considered), in this case, the transition probability can 
be written as, 



Wad-design) = y^[w a (l - n,)(A +A 1 (n i -i + n i+1 ) 

i 

+A 2 n i - 1 n l+1 ) + w d ni(D + Lq^-i + n i+ i) 

+D 2 ni-ini + i)]S(rii, 1 - n i )\lj^ i 8{n'j,n.j) (7) 
Here adsorption into site i occurs if rii = initially, 

with a rate controlled by prospective neighbors if A; ^ 0. 

The Kronecker delta for sites j ^ i excludes multiple 

transitions. 

The motion equation for coverage can be obtained by 
multiplying the master equation by occupation number 
rii and summing overall sites. Thus one can obtain the 
following expression for the coverage [1, H, [l(| [H, [H, HH 



— = w a [A E + 2A 1 NE + A 2 NEN] 
dt 

-w d [D Q N + 2D 1 NN + D 2 NNN] (8) 

With the same procedure the rate equations for the 
next correlation functions for three independent sites can 
be written as, 



= 2w a [(Ao + At)NE + {A 1 + A 2 )NEN] - 2w d [{D + D^NN + (D 1 + D 2 )NNN] (9) 

= w a [2A NEE + 2A 1 NEEN - (A + 2A X + A 2 )NEN] 

-w d [2D NEN + 2D 1 NNEN + (D + 2D X + D 2 )NNN] (10) 
= w a [2(Aq + A ± )NNE + (A + 2A 1 + A 2 )NEN + 2{A X + A 2 )NENN] 

-w d [(3L>o + 4L>i + D 2 )NNN + 2(D 1 + D 2 )NNNN] . (11) 

I 



w a A = w d D e 



(12) 
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w a (A + At) = w d (D a + D 1 )e-^ +v ^ 



and 



for i — 1,2, where 7 is a proportionality coefficient. 
Note that this linear relation was previously proposed in 
(13) references [HI, EH HI] , where A$ = 1 has been considered. 
The principle of detailed balance imposes, 



w a {A + 1A X + A 2 ) = w d (D Q + 2D X + D 2 )e 



-0(E„+2V) 

(14) 



w a and Wd cannot be fixed by detailed balance, be- 
cause they contain the information about the energy ex- 
change with the solid in the adsorption and desorption 
processes, which is not in the static lattice gas Hamilto- 
nian [22j. However, if one consider that w a = Wd — wo, 
and by comparison with the phenomenological expression 
for adsorption, one can identify it as 



Wq = Sq(T) 



PXa, 



(15) 



Sq(T) is the temperature-dependent sticking coeffi- 
cient at zero coverage and a s is the area of a surface 
unit cell. It contains the dynamic information about the 
energy transfer from the adsorbing particle to the solid 
which gives rise to its temperature dependence, for in- 
stance, an exponential Boltzmann factor for activated 
adsorption. It can be calculated only on the basis of 
a dynamic theory that accounts for the coupling of the 
ad-particles to the vibrational and electronic degrees of 
freedom of the substrate, and must be postulated ad hoc 
within the context of the kinetic lattice gas model. 

The long-range interaction introduces similar con- 
straints on the other coefficients. Each such constraints 
introduces two new coefficients. However, detailed bal- 
ance provides only half the number of relations to fix 
these unknown coefficients in the transition probabilities. 
Again, the static (lattice gas) Hamiltonian cannot com- 
pletely dictate the kind of kinetics pos sible in the system. 
As it is pointed out in references jl5l . [T9L |22| , any func- 
tional relation between the A— and D— coefficients must 
be postulated ad hoc, or calculated from a microscopic 
Hamiltonian that accounts for coupling of the adsorbate 
to the lattice or electronic degrees of freedom of the sub- 
strate. 

In the next section two different relations among the 
adsorption and desorption coefficients, are analyzed. 



Di=D 



■ 7e 



pv 



(17) 



and 



D 2 = D 



1 



7e 



20V 



2D X 



(18) 



Other possible choice of the parameters A; and Di, 
which fulfill the detailed balance and leads to a new for- 
mulation of the adsorption-desorption kinetics, is (case 
K2), 



and 



A 2 = jAi 



D 2 = 1 D 1 



(19) 



(20) 



Considering only nearest neighbor interactions one can 
replace (19) and (20) in the equations (13) and (14) and 
find the following expressions for the coefficients: 



A x = A 



(2 + 7)e-^-(l+7)e- 2 ^ y -l 
(2 + 7 )(l-e-' 31/ ) 



(21) 



and 



D 1 =D Q 



(1 + 7)e^ - (2 + 7) + ; 
(2 + 7 )(i_ e -/3V) 



(22) 



When adsorption and desorption coefficients are de- 
termined, the properties of the system such as adsorp- 
tion isotherms, sticking coefficient and thermal desorp- 
tion spectra can be obtained. 

However, to analyze correctly the consequences of the 
linear relations assumed above, let us introduce the fol- 
lowing identities 



III. LINEAR RELATIONS BETWEEN 
ADSORPTION AND DESORPTION 
COEFFICIENTS 

A linear relationship between the Aj and Di coefficients 
is proposed (case Kl) as: 



A 



A 
V 



(16) 



and 



E = NEN + 2NEE + EEE 



NE = NEE + NEN 



(23) 



(24) 



In such way, one can easily obtain alternative expres- 
sions for the equations (8)-(ll) as, 
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d£ 
dt 

dNN 
dt 

dNEN 
Jt 

dNNN 
dt 



w [{A )EEE + 2(A + A\)NEE + (A + 2A X + A 2 )NEN] 

-w [(D )ENE + 2(D + DJNNE + (A + 2D 1 + D 2 )NNN] (25) 
2w [(Aq + Ai)NEE + (A Q + 2A 1 + A 2 )NEN] 

~2w a [(A + D^NNE + (D a + 2D X + D 2 )NNN] (26) 
w [2(A Q )NEEE + 2(A + A^NEEN - (A + 2A 1 + A 2 )NENN - (A + 2A 1 + A 2 )NENE] 
-w [2{Do)NENE + 2(D Q + D^NENN - (D + 2D X + D 2 )NNNE - (D + 2D 1 + D 2 )NNNN] 



w [2(Ao + A X )NNEE + 3(A + 2A 1 + A 2 )NENN + (A + 2A 1 + A 2 )NENE] 
-wo [(2(£>o + A) + (D + 2D 1 + D 2 ))NNNE + 3(A) + 2A + D 2 )NNNN] 

I 



(27) 



(28) 



Note that EEE, NEE, NEN, ENE, NNE and 
NNN, as well as NEEE, NEEN, NENN, NENE, 
NNNE, NNNN and NNEE are mutually exclusive 
conditional probabilities, therefore each of the parenthe- 
ses in cqs. (25)- (28) must be positive. After some algebra 
one can obtain the following inequalities 



Ql ds = Ao>Q, 



Qlds = (A + A 1 )> 0, 



(29) 



(30) 



and 



Qdes = (Do + 2A + A) > 0. 



(34) 



y 
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Qlds = {M + 2A, + A 2 ) > 0, 


(31) 
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Q° des = A > 0, 


(32) 
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It it easy to see that eqs. (29)- (34) impose new restric- 
tions on the Ai and A coefficients. If these additional 
restrictions are applied to the two proposed kinetics, cer- 
tain values of parameter 7 are not allowed. This forbid- 
den values depend on the lateral interaction V, and this 
is shown in the diagrams of the Figs, la) y lb) for the 
Kl y K2 kinetics, respectively. 

Region I corresponds to those values of 7 and V which 
make positive all the coefficients in eqs. (29)-(34). Re- 
gion II corresponds to those values of parameters which 



FIG. 1: Diagram 7 versus V/UbT corresponding to a) Kl and 
b) K2 proposed kinetics. 



make Q\ ds > 0, Q des > 0, Q 2 ads < and Q des < (in- 
termediate region). Region III corresponds to those val- 
ues of 7 and V which make all the coefficients negative, 
Qlds < 0, Qlds < 0, Q des < and Q des < (forbidden 
region) . 

To calculate the exact solution for adsorption isotherm 
the two first rate equations (eqs. (25) and (26)) must be 
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FIG. 2: Adsorption isotherms for different values of the pa- FIG. 3: Sticking coefficients for different values of the param- 
rameter 7 and for two linear proposed kinetics, a) Attractive eter 7 and for two linear proposed kinetics, a) Attractive and 



and b) Repulsive lateral interaction. 



b) Repulsive lateral interaction. 



set equal to zero [13, EH . 

To calculate the sticking coefficient one used the fol- 
lowing definition 

S(6,T)=A E + 2AiNE + A 2 NEN (35) 

The exact solution for immobile TPD spectra is ob- 
tained by solving the rate equations for the first four 
correlations ((25)-(28)), where the adsorption terms are 
neglected [23, [5l|, [53| ■ In all TPD spectra the thermal 
desorption is considered as an activated process. Where 
the activation energy is 10 kcal/mol, the pre-exponential 
factor is 10 13 s _1 and the initial coverage is #0 = 0.9. 

The adsorption isotherms, sticking coefficients and 
TPD spectra for immobile adsorbate are showing in Figs. 
2, 3 and 4, respectively. 

For those values of the parameter 7 out of the allowed 
region, the first derivative of the isotherms presents a 
discontinuity as a function of the chemical potential, and 
both the sticking coefficient as TPD spectra take negative 
values. These behaviors are certainly anomalous. 

For values of 7 belonging to the allowed region the 
behavior of the observable do not present any inconsis- 
tencies. Although the physical meaning of the relations 
given in eqs. (16), (19) and (20) is not clear, it seem to 



be perfectly valid in view that they fulfill the detailed 
balance. However, it must be certain coherence in the 
behavior of the observables according to the lateral in- 
teractions. For example, for repulsive lateral interaction, 
the adsorption isotherms for monomers present a charac- 
teristic plateau at coverage 9 = 1/2 for enough low tem- 
perature. This is due to an ordering of the adsorbed par- 
ticles, namely, there is an alternation between particles 
and vacancies. For the same interaction the sticking coef- 
ficient will be below the straight line S(9) = 1 — 6, which 
corresponds to the Langmuir kinetics (null interaction) 
jl6j |. Finally, the immobile TPD spectra for high initial 
coverages will have three peaks, as is expected 0, [5o| . 
This behavior can not be reproduced by using both pro- 
posals, Kl and K2. 

Moreover, in the framework of proposal Kl, the three 
peaks in the immobile TPD spectra can be observed only 
in the limit 7 — > 0. However, this value of 7 corresponds 
to the Langmuir kinetics. 

Clearly, the principle of detailed balance is not enough 
to guarantee the correct behavior of the kinetics and 
deeper analysis must be done to choose the functional 
relation between the A— and D— coefficients. 
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FIG. 4: Immobile TPD spectra for different values of the pa- 
rameter 7 and for two linear proposed kinetics, a) Attractive 
and b) Repulsive lateral interaction. 



IV. INVERSE RELATION BETWEEN THE 
ADSORPTION AND DESORPTION 
COEFFICIENTS 

As is discussed in the last section, the relation among 
the adsorption and desorption coefficients must be care- 
fully chosen. In particular, linear relations between 
them, as is proposed by Kreuzer and co-worker in ref. 
13 HIJ are not valid in general and some values of 
the parameter 7 and V are not allowed. In fact, the 
observables obtained for those values of 7 and V are ill- 
behaved. In order to solve this kind of inconsistence, let 
us introduced the following relations between the A— and 
D— coefficients (K3), 



1 

An = , 

A) 



An + Ai 



1 



Dn+D! 



and 



A + 2Ai 



1 



D + 2Dx + D 2 



(36) 



(37) 



(38) 



Using the detailed balance expressed in eqs. (12)-(14) 
and after some algebra, one can obtain the following ex- 
pression for the A— and D— coefficients: 



A = exp[-f3E s /2] (39) 



D = exp[/3£ s /2] (40) 



A\ = exp[-/3E a /2] (exp[-/3V/2] - 1) (41) 



D 1 = exp[/?£ s /2](exp[/3V/2] - 1) (42) 



A 2 = exp[-/3E s /2] {exp[-0V/2] - l) 2 (43) 

and 



D 2 = exp[/3E s /2}(exp[/3V/2] - l) 2 . (44) 

Then, one can write the adsorption-desorption coeffi- 
cients as follow: 



and 



Q i ads =exp[-p(E s /2 + iV/2)], 



r des ^exp[f3(E s /2 + l V/2)}, 



(45) 



(46) 



with i = 0, 1, 2. These expressions for the coefficients 
allow an adsorption-desorption kinetics without physical 
inconsistencies. Particularly, in the framework of this 
formulation, Langmuir kinetics arises only for V = 0. 
In Figure 5 a), b) and c) it is shown the adsorption 
isotherms, sticking coefficients and TPD spectra for dif- 
ferent lateral interaction using the above formulation. As 
is observed, there are coherence between the behavior of 
the observables. 



V. SUMMARY AND CONCLUSIONS 

In this work, the adsorption-desorption kinetics in the 
framework of the one-dimensional lattice gas model with 
nearest neighbor lateral interaction is considered. The 
master equation approach has been used to derive the 
rate equations for coverage and higher correlations, which 
gives the equilibrium and non-equilibrium properties of 
the system. In order to obtain those equations, tran- 
sition probabilities are written in terms of the occupa- 
tion configurations. In such way, the principle of details 



8 




-6 -4-2 2 



10 12 14 



1,4 
1,2 

— \ 

o 

^0,8 
P 

<D 0,6 
0,4 
0,2 









\ 

\ 


V/k B T = -4 

V/k B T = 

- - -V/k B T = 4 




N ^^^^ 
X 





e 



.tS 0,08 

a 

=3 



-S 

CS 0,04 

H 
-a 

T3 



v = 


-2 kcal/mol 


v = 





- - -v = 


3 kcal/mol 




90 100 110 120 130 140 150 160 170 180 190 



T(K) 



FIG. 5: a) Adsorption isotherms, b) sticking coefficients 
and c) immobile TPD spectra, for attractive (dot), repulsive 
(dash) and null (line) interaction, using the inverse relation 
L3 



these parameters are wrong, or at least, the equilibrium 
and non-equilibrium observables present anomalous be- 
havior. A diagrams of the allowed and forbidden values 
as a function of the parameters V and 7, as well as the ad- 
sorption isotherms, sticking coefficient and TPD spectra 
for two different linear relations are showed. Three well- 
defined regions emerge from these diagrams in both cases: 
(i) The allowed region, where the adsorption-desorption 
kinetics is well behaved, i.e., the isotherms, sticking co- 
efficients and the TPD spectra are well behaved. It is 
also observed that Langmuir (Ai = 0) and interaction 
kinetics (Ai = —Di) belong to this region, (ii) The in- 
termediate region, where the first derivative of the ad- 
sorption isotherms are discontinuous. This results is not 
possible in a one-dimensional lattice gas with nearest- 
neighbor lateral interaction, (iii) The forbidden region, 
where the adsorption isotherms behave like in region II 
and sticking coefficients are negative. In order to solve 
such inconsistences, an inverse relation among adsorp- 
tion and desorption coefficients is introduced. In such 
way the equilibrium and non equilibrium observables do 
not present any anomalous or inconsistent behavior. 

As a general conclusion, the principle of detailed bal- 
ance is not enough to guarantee the correct behavior of 
the kinetics. The linear relations among the adsorption 
and desorption coefficients are not general and some val- 
ues of the lateral interaction V and 7 are not allowed. 
Inverse relations among the adsorption and desorption 
coefficients seems to be a better choice to determine the 
rate equations. However, deeper analysis must be done 
in order to choose a more general functional relation be- 
tween the A— and D— coefficients. 

On the other hand, both linear relations as well the 
inverse one belong to the so called soft dynamics, where 
transition probabilities factorize into a part due only to 
the change in the field energy and a part due only to the 
change in the interaction energy. 

Finally, it is important to emphasize that, even when 
the treatment has been done in one-dimensional system, 
the results seems to be general and they do not depend 
on the dimensionality. However, the extension to higher 
dimensions, the incorporation of next-nearest neighbor 
interactions, as well as of diffusion terms should be done 
in next future. 



balance imposes a set of restrictions on the adsorption, 
Ai and desorption, Di coefficients. However it deter- 
mines half of them. To overcome such difficulty differ- 
ent functional relations among those coefficients can be 
postulated, the simplest one is a linear relation. How- 
ever, it is observed that, when rate equations are written 
in terms of mutually exclusive conditional probabilities, 
additional constrains appear on Ai and Di coefficients. 
These are rather restrictive because not all the values of 
the lateral interactions, V and linear parameter, 7 are al- 
lowed. Moreover, the kinetics obtained for some values of 
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